December, 1993 



UBC-27/93 
ITEP-M5/93 

Evaluation of Observables in the Gaussian 
8 ■ N = oo Kazakov-Migdal Model 



M.I. Dobroliubov a 'Q' 2 , A. Morozov 6 , G.W. Semenoff °< 

N. Weiss a 'F| 

in ■ u 

a Department of Physics, University of British Columbia 
Cf) ' Vancouver, British Columbia, V6T 1Z1 Canada 

Q\ ' 

b Institute of Theoretical and Experimental Physics 
117259 Moscow, Russia 



X 



2 



ABSTRACT 

We examine the properties of observables in the Kazakov-Migdal model. We 
present explicit formulae for the leading asymptotics of adjoint Wilson loops 
as well as some other observables for the model with a Gaussian potential. 
We discuss the phase transiton in the large limit of the d = 1 model. One 
of appendices is devoted to discussion of the N = oo Itzykson-Zuber integrals 
for arbitrary eigenvalue densities. 
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1 Introduction 



The Kazakov-Migdal model (KMM) was introduced in |T| as an example of 
a gauge-invariant matrix model on a D-dimensional lattice. It consists of 
an N x N Hermitean matrix scalar field <&(x) which which lives on sites, x, 
and transforms in the adjoint representation of SU (N) and a unitary N x N 
matrix gauge field (link variable) U(xy) which lives on links, (xy). The 
partition function is 

Z= j d${dU] exp | -N^2TtV{®(x)} + N £ Tr$(x)U(xy)$(y)rf(xy)\ 
\ x <x,y> j 

(1) 

where V[<&] is a self- interaction potential for the scalar field. Due to the 
absense of a kinetic term (such as, for example, a Wilson term) for U, one 
can integrate out the gauge variables in ([!]) exactly.^ The resulting effective 
scalar theory contains iV degrees of freedom, the eigenvalues of the matrix 
<&(x), whose quantum fluctuations are suppressed in the large N limit. This 
makes it possible to employ large N techniques. The field <&(x) plays the 
role of a master field and the density of its eigenvalues, p(4>), can in principle 
be found by solving a saddle-point (master field) equation. In this sense the 
model is regarded as being solvable in the large N limit. However, explicit 
solution of the master field equation remains one of the major technical 
obstacles in this program. 

With the assumption of a translation invariant master field, Migdal @ 
has reduced this problem to finding the solution of a singular highly non- 
linear integral equation. Some further progress has also been made using 
the method of loop equations [|J. In both approaches, technical problems 
prevent one from obtaining very many analytic results. For instance, explicit 
solutions for the spectral density of the scalar field are thus far known for only 
two kinds of potentials: the Gaussian one (V[$] = ) solved by Gross 

U and a potential related to the Penner model, which was investigated by 
Makeenko ||. Furthermore, even for a Gaussian potential where the spectral 
density is known explicitly, only expectation values of products of scalar fields 
on a site and one-link correlators of the gauge field || have been computed 

1 The exact integration of the gauge fields from (TJ) uses the Itzykson-Zuber formula Q 
and can in some sense be regarded as a way to avoid solving the long-standing problem of 
summation of planar diagrams in lattice gauge theory. 
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so far. In this Paper, we shall use these previous results to present some 
explicit formulae for observables in the Gaussian KMM. In particular, we 
give a caluclation of the expectation value of the adjoint Wilson loop as well 
as the gauge invariant two-point correlator of the scalar field. We regard this 
as an intermediate step in understanding the continuum limit of the KMM 
in the general case. 

It was originally hoped that the Kazakov-Migdal model would possess a 
continuum limit which would be a gauge theory, such as QCD. It was also 
clear from the outset that this could not be realized in a straightforward way. 
Because the action of the model ([I]) does not contain the kinetic term for the 
gauge field (this was exactly the fact which allowed one to integrate out the 
gauge fields), the KMM posesses an additional local (gauge) invariance 
[0] which implies that all of the usual Wilson loops vanish unless they have 
vanishing area. This can be interpreted as an infinite string tension and 
results in "ultraconfinement" , where no propagation of color is allowed at 
any distance scales. An ordinary continuum gauge theory does not have 
such a property. 

Thus, if the KMM is to describe continuum gauge theories, then either 
there must exist a more sophistocated continuum limit M, or the model 
itself should be modified in a way which would break the local Zn symmetry 
explicitly (yet preserving its solvability) ||, [iOfl . 

Irrespectively of its relevance for the description of D > 2 Yang-Mills the- 
ories, the KMM has proven to be very interesting from a purely theoretical 
point of view as an intermediate step between generic matrix models which 
involve full unitary matrix integrals, and the well studied "c < 1" matrix 
models in which the unitary matrices ("angular variables") completely de- 
couple The KMM is furthermore directly relevant to a certain sector of 



c = 1 string theory where it indeed has a phase transition and one can 
discuss its continuum limit. 

The choice of physical observables in the KMM is not entirely obvious. 
Due to the existence of the additional Z^ gauge invariance, all of the con- 
ventional Wilson loops vanish. Note, that in the large N limit and in the 
mean field approximation, the unitary integrals in the KMM can be reduced 
to single link integrals of the type 

C(n|<&, 9) oc / [dU]U nn . . . U injn U] h ...UU e ™ w ** , (2) 
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They must involve the same number of U's and W's in order to be invariant 
under the local Zn transformation. The simplest gauge invariant observables 
which can be constructed out of the C's are the adjoint Wilson loop and 
the filled Wilson loop. Both of these involve only the averages of \Uij\ 2 i.e. 
C(l|$,$). (As we mentioned above, even for a soluble model such as the 
Gaussian KMM C(n\<&, 5>) have been evaluated explicitly only for n = and 
n = 1). 

Before presenting these formulae we make some general comments about 
the integral (§) for arbitrary n. These integrals fall into the general frame- 
work of applicability of the Duistermaat-Heckmann theorem, and, in prin- 
ciple, they can be explicitly evaluated (see || and references therein for 
details). However, existing explicit expressions for n > 1 [13] are too com- 



plicated to be used effectively when N — > oo, which is the limit of interest 
in the KMM. What is required is an expression for C(n|$, $) in terms of the 
densities p(0) = 4 Y,iLi ${<t> — 0i) an d p{<f>) of eigenvalues of the matrices $ 
and $. This would allow taking a smooth large N limit. In this limit, dis- 
crete sums over the indices become continuous integrals over the eigenvalues, 
J2i — * I d(pp((f>). Expressing all quantities in terms of densities implies also 
relabelling of tensor structures like C(n) for n > 1. Indices should 
be substituted by 0j, <pj . . . according to the rule Cfc ...«^ im = N m ^ 2 Ci 1 ,,,i m , so 
that for example 

N x 

^QjCjk = -j^zZC^-C^ = / d^p{4>)C^C^ k , (3) 
i=i *i 

At the moment there exist only a few methods for the evaluation of C(n\p, p) 
with n = |nj and n = 1 | (for details see Appendix B). Though they 



have the potential to work in a general case, these methods have not proven 
to be very efficient in practice. Until now the quantities C(n\p,p) have been 
computed only for two cases: for n — 0, 1 for the semi-circle |5|, |J and n = 1 
for the "Penner-like" |J distributions of the eigenvalues. (The latter is more 
complicated and will not be considered here.) 

With the assumption that the mean field is spatially homogenous, the 
scalar field in the Gaussian KMM has the semi-circle spectral density 



P (0)# = p(0)# = P ^)d4 = lp- M < -j=) ( 4 ) 
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For this spectral density the explicit results for the Itzykson-Zuber integral 
C(0\pfj,, p^) || and correlators C(l\p^,p^) f| are as follows 



C(p)= lim -llog / [dU]e 



Tr7V<t>[/<I>[/t 



C{p) = ^ TT ~ A -^Uo g ^ TT ~ A + » ; (5) 
vw 2/i 2 & 2/i v 7 



C^fr) = lim NCM = \ + ^f+^ . (6) 

In practical calculations it is more convenient to use the "normalized" semi- 
circle distribution, p s ((f)), 4> = 3p,4>i which has support [—1,1] and 

p s (0) = Vl-0 2 , p(0)# = 2p s (0)#, (7) 

7T 

so that 

„, . coshr — sinhr 1, cosh r + sinh r 
^( r ) = 7T—L n lo § 



2 sinh r 2 2 sinh r 

sinh r(cosh r + sinh r) 

d 2 + /3 2 — 2d/3 cosh r + sinh^ r 



- 2C a =& & ,P=&p- ~o , no 7^~~n i , • i 2 ' ( 8 ) 



where sinh r = — , cosh r = — — — . 

2' 2 

(sometimes we shall use notation C & &(t) instead of C & ^(p), which should not 
cause any confusion). 



2 Observables in the KMM 

The set of observables in the KMM is restricted by the local Z7v-invariance, 
which, for example, forces averages of all the fundamental-representation 
Wilson loops to vanish and results in ultraconfiment. Non- vanishing observ- 
ables contain an equal number of U and W matrices at every link,0 examples 

2 Since we consider the large- N limit, we neglect "baryonic" observables, containing N 
or more U -matrices per link without any 
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are shown in Fig.l. Dots can appear at every site and denote insertion of 
any number of $ or $ operators. One can account for all possible insertions 
by using the generating function , at every site. Occasionally some 
double-lines can become quadruple etc at some links, then the knowledge of 
higher-order correlators C(n) with n > 2 is required for evaluation of aver- 
ages. We restrict ourselves to consideration of the non-degenerate situations, 
where double lines never overlap. In such cases the answer is further sim- 
plified || by the occurence of 6*6 ^ in eq.(|). It allows us to consider every 
observable as associated with a graph T which ascribe a single index i x to 
every site (i.e. effectively substituting double lines by single lines, with CV,- 
playing the role of the propagator, and (X x — (fi^)' 1 - that of the vertex): 

^ iy -n r+# (, e r)- # ( <a;2/ > e r) jj j P^tk JJ C My , (9) 

where, as usual, nr is the number of traces, contained in the definition of 
Or (for discussion of different normalization prescriptions see ref.||). Note 
that, with our normalization of C and p, each site contributes a factor of N 
and each link a factor of 1/N. It will turn out that the product of integrals 
in the right-hand-side of (9) is well-defined and iV-independent. With this 
normalization of operators, the only ones which are non-zero in the infinite N 
limit are tree-like configurations - these always have one index sum (rip = 1) 
and also #(x G T) — #((xy) G T) = 1. Every loop in a configuration will 
only increase ${x G V) — ^((xy) G T), making it greater than one. 

At least in principle, the integral on the r.h.s. can be evaluated for any 
graph, provided C<j> x <t> y is known for the given p. At the moment, however, an 
explicit expression is available only for p given by the semi-circle distribution. 
In D = 1 there are only two types of allowed graphs: segments and circles of 
length L which, as we shall show, can be relatively easily computed for the 
semi-circle distributionP]. This is the subject of the next section. In higher 
dimensions the full set of observables contains more objects, like filled Wilson 
loops 0, which are not considered in the present paper. 

3 For lattices which form "rare nets" , so that the theory remains essentially 1- 
dimensional, more sophisticated graphs are also of interest. We, however, ignore this 
possibility. 
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3 Evaluation of Observables 



3.1 The Adjoint Wilson Loop 

The Wilson Loop is obtained by considering the trace of the path-ordered 
product of the link operators for the links which occupy a contour T on the 
lattice 

W[T) =tr [] U ( x v) ( 10 ) 
<xy>er 

Because of the symmetry, the expectation value of this operator vanishes 
unless T has vanishing area. An example of an operator which is invariant 
and also a physical interpretation is the adjoint Wilson loop, 

^[r] = ^(|^[r]| 2 -i) (ii) 

The appearance of the constant term in this definition can be understood if 
one recalls^ that the product of two fundamental representation link opera- 
tors contains both the adjoint and scalar representations, the latter of which 
must be subtracted: {UijU^A = U^U^ — SuSjk- 

In the mean field approximation of the KMM, the adjoint Wilson loop 
can be evaluated using the one-link expectation values of the the products 
On = (W), 

M = i 5 (tr(^)-l) = i 5 (i:cf n -l) (12) 

where q are the eigenvalues of the matrix Cy and L\T] is the length of the 
contour T. 

From their definition it is easy to see that the matrices CV,- are real, have 
positive entries and, when the master field is homogeneous, they are also 
symmetric. Furthermore, they obey the sum rule 

N 

E^ = 1 ( 13 ) 

i=l 

This sum rule implies that, for any N, Cij always has one eigenvalue which 
is one, with eigenvector (1, 1,1,..., 1). 

4 We thank Yuri Makeenko for a discussion of this point. 
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It is also easy to see that all other eigenvalues lie in the interval [—1,1]. 
Indeed, let us consider the set of real matrices with positive entries, satisfying 
the sum rule . Obviously, all matrices from this set have their traces 
bounded from above by N. Further, this set is closed under the matrix 
multiplication: if matrices A and B belong to this set, then so does the 
matrix AB. Take now any matrix C from this set and consider its power 
C n with n even. If there are eigenvalues of the matrix C, absolute value of 
which is greater than 1, then for big enough n one can make the trace of C 
arbitrary large, thus violating the mentioned upper bound. Q.E.D. 

Let us recall that a rectangular adjoint Wilson loopp (whose lengthes in 
the space and time directions are L and T, respectively) in the limit T — > oo 
can be interpreted as the energy of a pair of mesons with separation L, 

E[L] = hm i In (n 2 < W A [T] >) (14) 

Once 1 is subtracted in the sum on the right-hand-side of (|12j), the energy is 
dominated by the largest remaining eigenvalue of C. 

In the following we shall compute W A in the large N limit for the Gaussian 
Kazakov-Migdal model. First, however, let us consider the case of iV = 2 
which was previously analyzed in ||. There, even though the mean field 
approximation for the scalar field is uncontrolled (since N = 2 is not large), 
it has been argued [15| that mean field theory can give accurate results. If 



the mean field eigenvalues of the 2x2 scalar matrix field (which can be 
taken as traceless) in this case are and —<f), it was shown in || that 

„ _ 1 - 1/40 2 + e- 4 ^/40 2 

Wl — ^22 — 



1 - e- 4 ^ 2 

r r V40 2 -e-^ 2 (l + l/40 2 ) 

O12 - C 2 i - _ e _^ 2 [Lb) 

The eigenvalues are 1 and coth(20 2 ) — 1/20 2 . Note that the second eigenvalue 
varies between zero and one as <j) 2 goes from zero to infinity. The energy of 
the meson pair is 

E = -2 In (coth(20 2 ) - 1 /2<p 2 ) (16) 

which is positive and goes to zero at the "critical" value of <p 2 — > oo. Note 
that, in this approximation the mesons have no interaction energy, i.e. (|I6|) 
can be interpreted as twice the meson mass. 
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Now let us proceed to the large N limit in the case of a Gaussian potential. 
We start with the basic relation, proven for the semi-circle distribution in 
Appendix A: 

C^{T X )C h {T 2 )p s 0)dp = C^(n + r 2 ). (17) 
Therefore for the adjoint Wilson loop of length L we find 

(WaF]) = (j daip(ai) . . . da L p(a L )C aia2 (r) . ..C aL _ iaL (r) 
= Jp(f d"ip("i) • • •d6cLp{pt L )C &l&2 {T) . ..C &l _ 1&l {t) - 1 



1 

iV2 



dap(a)C 'aa(Lr) 

1 e- Lr 
iV2 1 - e~ LT 



Note that, as noted in the previous section, the normalization of the adjoint 
Wilson loop is such that it vanishes in the large N limit. We can regard (IT3) 
as the leading asymptotics. 

The interaction potential for a pair of mesons is obtained from the free 
energy, 

E(L) = 2r (19) 

which one can interpret as just the sum of the meson masses (each equal to 
t) without any interaction between mesons. 

There is an apparent phase transition at r = 0. Recall that ||, for the 
Gaussian potential V[<p] = ^4> 2 , so that the "bare mass" of <p is m2 ~ 21?, 
and r is given by 



m \D-l)+D^m±-A(2D-l) 

sinnr = ; — ; 20 

2(2/?- 1) v ; 

or, equivalent ly, 

m 2 = 2Z?e~ T + 2sinhr^2Z? + 2(l-Z?)r (r ~ 0) (21) 

As has been observed by Gross ||, the phase transition can be approached 
through the physical region, m 2 > 2D, only when D < 1. 
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3.2 Other Observables 

It is straightforward to compute the expectation values of a variety of other 
observables. Here, we shall present a few examples. 

It is very easy to find the expectation value of powers of the scalar fiels 
$ on a single site, using the generating function 

E x = (tr ) = V il^LZ . (22) 

This generating function for the semi-circle spectral density was found by 
Gross || (in this case odd momenta are vanishing since the spectral density 
is even) 



& = /T37r = T-V— -"=A[ 1+ g-(FM)r^) I • (23) 

so that (tr$ 2 ) = l//i and so on. 

For an adjoint loop of the length L with the insertion of the operator — j| 
into one of the products II U at the single site we obtain: 

If 1 

(0(oL)) = — / donp{oti) . . .da L p{a L )- C aia2 (r) . . . C aLai (r) ■ 

l Vp r l 



/ dap(a)- -Ca&{Lr) 

J A — n 



iV 2 2 J ry '\-a 
1 e- Lr 1 I 



(X--(e LT + l)E x ) (24) 



N 2 l-e~ LT X 2 - Jcosh 2 (Lr/2) v " /i' 

where is defined in the Appendix. The coefficient in front of A -1 in the 
large-A expansion of (|24|) reproduces (|l8l). 

Similarly, for a segment with ^-insertions at its ends: 



f 1 1 

/ daxpioii) . . . da L p(a L )- C ai 

J A — ct\V — a L 



!\ 2 \ T ) ■ ■ ■ ^ot L - 1 a L \ r ) 



= t j i& ^ 43 m c (Lt) 

A — a v — p 

XE V + uE x + 2 sinh(Lr) + (cosh(Lr) + sinh(Lr))( 2sinh(Lr) E^ A - XE X - vE v 



z/ 2 + A 2 - 2 cosh(Lr)z/A + J sinh 2 (Lr) 
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It is worth noting that, for real r, the expressions ( |2~3] ) and ( p5|) are never 
singular for any A and u, as they should. 

Also, note that this quantity is of order one, rather than 1/N 2 . This is a 
result of the fact that this observable is tree-like, whereas the previous two, 
which we considered, had loops (cf. the discussion after the equation (Q)). 
The gauge invariant <fi — <fi correlator, which is a flux tube with the field $ at 
the edges, can be extracted from this result by taking the leading, order of 
1/A 2 , l/z> 2 , asymptotics of (p5|) as 

(tr <p\{U(i) \{U ] ) = -e~ LT (26) 
r — r ft 

We see, that it has very simple form, with r being the correlation length 
and agrees qualitatively with our interpretation of the result (|T9|) . When the 
length L equals zero, so that the flux tube vanishes shrinking to a point, the 
expression (|26|) reproduces the result for the average of tr$ 2 . 

4 Conclusion 

The effective field theory for the eigenvalues of the scalar field $ is classical in 
the large N limit, and the classical configurations are obtained by solving the 
saddle point equation. It is a remarkable feature of the Gaussian Kazakov- 
Migdal model that, in addition to determining the classical scalar field, one 
can take the fluctuations of the gauge fields into account exactly. 

From the trivial counting of powers of iV we find that, with their conven- 
tional normalization, the only operators with non-zero expectation value in 
the limit N —>■ oo are those of tree-like configurations of tubes of glue with 
powers of the scalar fields inserted at the ends of the branches. 

While the expectation value of the Wilson loops with non-zero area is 
always zero in this model, the expectation value of the adjoint Wilson loop 
is proportional to 1/iV 2 and obeys the perimeter law, with the correlation 
length equal to 1/r, which is related to the bare scalar mass m 2 — 2D and 
the coordination number of the lattice 2D according to (21), (22). In the 
continuum limit (which is possible only in D = 1), this corresponds to non- 
interacting mesons with the mass equal to r. 

For a flux tube with scalar field at each end, the expectation value is 
proportional to the exponential of the length of the tube divided by the 
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correlation length, which is again given by 1/r. We also obtained the general 
formula for the expectation value of the flux tube with arbitrary powers of 
the scalar fields at each end. 



APPENDIX A 

This appendix contains some formulas, relevant to the proof of eq.([T7|) as 
well as for other calculations with semi-circle distributions. 

We repeat the definition of the "normalized" generating functional for 
the averages (Tr$ 2fc ) 

- _ f p,0)dp „ r— - i / (2k-iy\ ( i Y' 

* = J -jZj- = ~ ~Tk + TjT \W) 

Eiy = cosh 9 — sinh 9 , where 7 = cosh 9 . 

Let us rewrite the correlators (^) for the "normalized" spectral density in 
the form 



g sinh r(cosh r + sinh r) sinh r(cosh r + sinh r 

6/9 ~ a 2 + (3 2 -2aj3 cosh r + sinh 2 r " - a+)0 - a J) 

whereQ 



a± = a cosh r ± \J a 2 — 1 sinh r = cosh(r ± t) , a = cosht (28) 



From (p7|) and (|28|) we get 

£ a± = cosh(r ± t) - sinh(r ± t) (29) 



(here the choice of the signs corresponds to that of Eq. (|23|)) 

The first applicaiton of this convenient parametrization is to check ex- 
plicitly the normalization condition (Rf) (we denote a = cosh 9) 



J d(3p((3)C a/ 3 = 2 J d(3p((3)C & p = sinhr(coshr + sinhr 
sinh r(cosh r + sinh r 



a_ — a + 

cosh(r + 9) — sinh(r + 9) — cosh(r — 9) + sinh(r 



cosh(r — 9) — cosh(r + 



3 Note that since a 2 < 1, t is complex 
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Let us now turn to the convolution of two C's. 



sinh 7;(cosh T; — sinh 7$) sinh (cosh Tj — sinh 



~ P-i P+j ~ P-j 

{ (3+i — (3+j (3-i — /3_ j f3-i — f3 + j f3 + i — P-j J 
sinh Tj(cosh r« — sinh Tj) sinh Tj (cosh Tj — sinh Tj) ( ^ T\ + 7j + 0i + 2 

. , Tj + 7} , , Ti + Tj -6 1 + 6 2 Ti + Tj + 01 - 2 1 

+ coth - coth coth > = 

2 2 2 J 

sinh(ri + r 2 ) (cosh(ri + r 2 ) + sinh(ri + r 2 ) 
PI + M - 2p x p 2 cosh(Ti + r 2 ) + sinh 2 (ri + r 2 ) ' 

which proves Eq. (|T7]) . 



APPENDIX B 

DMS-Matytsin theory of the N = 00 Itzykson-Zuber integrals 

This appendix contains description of the currently available indirect 
methods to evaluate the (logarithm of the) Itzykson-Zuber integral C(0\p,p) 
and the first Itzykson-Zuber correlator C a g(l\p, p) for given eigenvalue den- 
sities p((p)d(p and p(0)d</>.Q 

As it often happens in the field theory, evaluation of the correlator C a p jp 



is somewhat simpler than that of the "free energy" C [H], and the methods, 
used in the two cases, though similar in many respects, are still complemen- 
tary rather than identical. 

Evaluation of C a/ 3 [|J 

Despite it produces the answers for individual (one-link) Itzykson-Zuber 
correlators the method of [ffl actually works within the framework of the 
KMM and makes use of the approach of the loop equations. 

6 In order to avoid confusion we emphasize that "bar" does not mean complex conjuga- 
tion, p(</>) and p{4>) are just independent functions. Complex conjugation will be denoted 
by "*" in what follows. 
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The main object of this approach is the following generating function [] 

W ^(7) (— ■ (30) 

It satisfies the following loop equation ||] 

(E a + (3)G a ? - E a / &L L W G »f> = L{a)Gap _ # ( a) . 

J Am a — UJ 

(E P + a)G aP - E/3 ± I fo L f) G °» = ~ mG _ £ M (31) 

J 2m p — uj 

where 

so that Disease, = 2mp(P), while Disc^E 1 ^ = 2mp((3). Functions L(u) and 
L(a;) are some functions. On the support of the corresponding spectral den- 
sity, these two functions are sums of the derivative of the logarithm of the 
Itzykson-Zuber integral and the real part of the generating function E(E). 
Both L and L are analytic in the vicinity of the supports of p and p and 
in principle can have cuts and singularities elsewhere. For the sake of sim- 
plicity we assume below that these singularities are located at uj = oo and 
uj = oo, and neither L nor L have cuts. These two functions are a very 
natural starting point for the DMS procedure. If one assumes an ansatz for 
L and L, then the search for the corresponding E a and C a p can be reduced 
to an algebraic problem [|], [J. Here, however, we will offer another variation 
of DMS approach, which is closer to the Matytsin's ideas of computing the 
Itzykson-Zuber integral and considers instead of L, L their combination with 



E a: E a (|33) 



The functions Rp{a) and R a ((3), enetering right hand sides of Eq. (|31|) , 
are defined as the contributions to the integrals (|3l|) from the residues at 
infinity and possible singularities of the functions L and L. According to its 
definition Rp{a) has no discontinuities as a function of a, Disc^-R^a) = 0. 



7 Note that the second scalar factor contains <I>. 
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Moreover, as a function of a it has its singularities only at those of L(a) 
(which are actualy singularities of gR(a) - see eq. (|33|) ). As a function of 
(3, Rp(ot) is analytic outside the support of p, where its imaginary part has 
a finite jump. Under our simplifying assumption this means that Rp(a) is 
a polynomial in positive powers of a with /3-dependent, everywhere finite, 
coefficients.^ 



From (|31|) it follows that 



Eg - Rp(a) _ Ep - R a (f3) , . 

P - g R {a) a - g L {(3) 



where 



9r(oc) = -E a + L(a) 
g L (j3) = -E f) + L(j3) . (33) 

According to their definitions E a ~ -, while G a p ~ — as a — > 00. 
Then C ai g is nothing but a double discontinuity of G a/ g: 

Disc a Disc^ Disc /3 i? /3 (a)/27rip(/?) 
Q/3 " 27rip(a) 27rzp(/3) 0/3 = ((3 - g R (a))((3 - g R (a)) ' 1 j 

where we have defined on the cut 

1 ~ 1 ~. 

g R (a) = L(a) - -V (a) + iirp(a) , -V (at) = Re£ a , (35) 

and "*" denotes complex conjugate. Alternative expression is 

Disc a fi a (/3)/27rzp(a) 

Q " " (« -<to(fl)(« -£(/*)) " 

Assume that the correlator C a p can be naively analytically continued (i.e. 
just by using instead of a and (3 two arbitrary complex numbers) from the 
cut. Since for any given a the positon of the poles of the expressions (|34]) 
and (l36|) should concide we find that 



g L (g R (a)) = a, g R (g* L ((3)) = (3 . (37) 



8 If L(a) has extra singularities at some point i/, Rp{a) will also contain a contribution, 
which is polinomial in (a — v)^ 1 ■ This correction is easy to be accounted for in our 
reasoning below. 
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Thus we recover the functional form of the Master Field Equation in the 

) . In order to do 



Maty t sin's approach . 

It remains to define the quantity in the numerator in 
this let us return to eq.(|32|) and rewrite it as f\ 



Discg-R^a) 
*™p{P) 



{(3-g R {a)) 



Disc^G, 



aj3 



P - 9r(<x) 



a 



7 



(3? 



Now it is time to recall that Rp{a) and thus Disc^-R^a) is a polynomial 
in positive powers of a (provided g R and L are singular only at infinity). This 

at the r.h.s. by the part of its asymptotics in 



allows to substitute 



P-9r(<x) 
a ~ 7 

a, containing only its positive powers 



P - 9r{&) 




9r(oc) 


a — 7 


+ 


a — 7 



ak l l(J kl 



(39) 



k,l>0 



i.e. by a finite polinomial in a and 7 with coefficients a^i = cw{<?.r}, totally 
defined by the shape of the function gn(a). 
Thus 



where 



DiscpRp(a 
2™P~(P) 



Mm 



k,l>0 



(40) 



yp( T )d T c^. (41) 

Note that it follows from the normalization condition for C a p that 

M (P) = 1 (42) 

Finally, from ([34]) we obtain: 

Ek,i>oa k Mi(P)a H {g R } 



(P - g R (a))(P - g* R (a)y 



(43) 



We use the fact, following from the definition (|30D, that 



1 

2nl 



Disc^Ga 



p{l)d-y 
a — 7 



C 7/9 p(/?). 
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and Mi((3) in this formula can be obtained from solution of a finite system of 
equations, which arises after ([43D is resubstituted into (|4lD with the condition 

This provides a complete solution for the problem of evaluation of C a p, 



provided one starts from any adequate (i.e. satisfying (j37D) pairs of functions 
9l,r[p)- We did not prove here that any such pair of functions gL,R.{ct) 
provides an answer: we rather proved that any answer has this form with 
some gL,R(ot)- It is clear, however, that further restrictions on the shape of gn 
should not be imposed: otherwise it would be impossible to suite arbitrary 
densities p and p. Therefore the alrorithm is obliged to work. 

Unfortunately, in all but the simplest cases the analytic structure of L{a) 
and/or E a is very complicated, which hinders any practical calculations. One 
can see it by following. Take any function g* L , which has negative imaginary 
part when its argument belongs to some real interval (support of p). Then 
consider the inverse of this function and call it gn (it must have an interval 
on the real axis, the support of p, where its imaginary part is positive). In 
most cases the analytic structure of gn turns out very involved. 

For illustrative purposes we turn to the bi-semi-circle example, where 
both p and p are semi-circle. We start from the pair of functions 



lp 2 + Ap/p U^a 2 
9l{ol) = - 1 a + \ — fj, 



4 



p 2 + 4/x/ p 



/ U 2 CV 2 

9r(P) = 1 2 a + V 4 ~ P (44) 

First let us find the coefficients a^i entering the power expansion 

1 



ooo = p 2 + Ap/p + p) , a k>0 ,i >0 = (45) 



The general result for the correlator (|43D then takes the form 

c nip p-) = - V/x 2 ^ 2 + 4pp + pp 

a/3 M ' ^ 2 pa 2 + p(3 2 — \/p' 2 p 2 + 4ppa/3 + pp ' 

which coincides with the result of |4[]. 

Knowledge of the correlator allows one in principle to compute the Itzykson- 
Zuber integral itself. To demonstrate that let us introduce an auxilary de- 
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pendence of C (n) on some parameter t in the following fashion 



Ct(0) = ^- 2 log j[dU}e m ™* u \ 

Then for its derivative we find 

dC t (0) J[dU}j ! Tr<5>U$We NtTr ' s ' mu ^ 



(47) 



dt (1) 

dapt{a)ad(3p t ((3)(5C af j(t) 



(48) 



As an example let us consider again the case of the bi-semi-circle distribution. 
The dependence on t can be absorbed into the redefinition \i and/or p,. It is 
convenient to consider the symmetric redefinition 



Then 



dC t (0) 
' dt 



2t 2 



/JIfI(y/ flfM + 4t 2 + y/JIJi) 



(49) 



(50) 



Integrating the last line one obtains the generalization fl5TD of the Gross' 
result 



c(o\ P „ Pll ) = - l - log ^±M 



(51) 



Matytsin's theory 



A more direct way of calculating the Itzykson-Zuber integral C(0), which 
does not require the knowledge of the correlators C(l), was developed re- 
cently in 



14]. It has been proved that 



C(0\p,p)=S(p,p) 
p{<t>)p{(f>') log(0 - <P 



p(<P)<P 2 d<p 

,_l 
2 



p{<t>)4> 2 d<p- 



p(<j>)p(<j>') log(0-0')W- 
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Integrals are over supports of p((f>) and p((f>) and 



S(p,p) = \j^J d<f>[lmf(4>,t)} ([Re/(0,t)] 2 + i[lm/(0,t)] 2 ) . (53) 
The function f(<f),t) is a solution to the "Hopf equation''^ 

Tt=% (54) 

with the initial condition 

/(0,t = O) = ^(0)-0. (55) 

It is then a consequence of the Hopf evolution that 

/(<M = l) = -S L (0) + 0. (56) 

The function f(4>,t) is in fact determined from a algebraic equation, because 
the differential Hopf equation (|54D is in fact explicitly integrable. Its generic 
solution is given in the parametric form: 

(f> = a + tF(a), 

f{<f>,t) = F{a), (57) 

with arbitrary function F(a). The shape of F(a) is dictated by initial con- 
dition (1551): 



t = 0)=g R (<p)-<p. (58) 
Given this function, one finds a(4>, t) from the algebraic equation 

<p = a((p,t) + tF (a((p,t)) , (59) 

10 It is a remarkable equation, the first one in the "quasiclassical KdV hierarchy" , which 
plays an important in the theory of topological Landau- Ginsburg models and Generalized 
Kontsevich model. In Itzykson-Zuber theory 1/t plays the role of the coefficient in front 
of the action, jTrQU&W . It is interesting to understand the role of other KdV times in 
this context. The "action" S(p, p) in ( |53| ) is not an action for Hopf equation - it is rather 
one of the (quasiclassical-) KdV Hamiltonians. We note in passing that the action for the 
Hopf equation is not unique: see ref. Jlt| (where it is refered to as Bateman equation) for 
discussion of the corresponding universality structure. 
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and then obtains 

f{<f>, t)=F («(0, t)) = g R («(0, t)) - t). (60) 

In order to make this (exhaustive) description a little bit more transparent 
we turn to the semi-circle case. 



p(a) = p(a) = Pfl (a) = ^-Jp- 



vV 2 + 4 . /— f <• cosn T 



g R (a) = a h inp^a) = y/Ji f a ^ + iy/l - a 2 



yTZ smhr smhr 



Then 



= t = 0)= g R {<j>) -<j> = y/ji ( C °^J hr \ + i VI " 2 ) • (62) 



Given such F (^) is actually a quadratic equation for a((f>, t), which is easily 
solved: 



0(1 + t(coshr - 1)) + iyi + 2t(l - t)(coshr - 1) - 2 
° ( *' t} = l + 2t(l-t)(coshr-l) (63) 



and substituting this into (60) we get: 



s9^T = i(i_2tU + iJl + 2*(l-t)(coshr-l)-^ 2 , s 
/(« = V^^ l + 2 f (l- f )(co shT -l) — ■ < 64) 



(One can now see that /(</>, t = 1) = ^//J l^- ™^ 1 ^ + zy 1 — 2 J is indeed 

equal to — #l(0) + 0, in accordance with flSEj).) 

Integral (p3|) for S(p fl ,p IJ ,) is now easy to evaluate and we get: 
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To find the Van der Monde determinant in the case of the bi-semi-circle 
distribution, we first note that it is easy to calculate the following integral 



p s , m (a)p SjAt2 (/3) log(>a - (3)dad(3 = 

(66) 



- 2 <S' Tf -|(l + log/i 2 ) 



z 2 >ia ia-lfi + ioga 



Thus, the logarithm of the Van der Monde determinant for the semi-circle 
distribution is 

pMPnW) lQ g(0 - 00 W = ~\ log/i - \ (67) 



Recalling that 



P Mrd<P = -, / PMfW = - , (68) 

we reproduce Gross's answer || for semi-circle distribution, 

C, = ^±i^_I log vZ±l±^. (69) 

2/i 2 2yU V ; 



We aslo note in passing that the item S(p, p) at the r.h.s. of (|52|) can be 
distinguished from the other terms for the following reason. Exact answer 
for the Itzykson-Zuber integral is: 

1 f r Tr<E>i/3>c/t detjje^ 1 ^ 



[dU)e 



V N J l J A(0)A(0) 

iTr$ 2 + iTrl> 2 

= _ A e t.. e ~l^~4>j) 2 C7Q) 

A(0)A(0) det ^ e (7Uj 

It is the last determinant at the r.h.s. which corresponds to e s( - p '^ in the 
large N limit. This determinant tends to unity for very broad distributions 
when the average distance between adjacent eigenvalues is large. In the 
example of semi-circle distributions this corresponds to p, p — > and indeed 
S(p,p) ~ 0(1), while the other terms are ~ 0(l/p,\ogp). 
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